Third order spectral branch points in Krein space related 
setups: VT— symmetric matrix toy model, 
MHD a 2 — dynamo, and extended Squire equation *) 

UWE GUNTHER AND FRANK STEFANI 

Research Center Rossendorf, Department of Magneto-hydrodynamics, 
P.O. Box 510119, D-01314 Dresden, Germany 

The spectra of self-adjoint operators in Krein spaces are known to possess real sec- 
tors as well as sectors of pair-wise complex conjugate eigenvalues. Transitions from one 
spectral sector to the other are a rather generic feature and they usually occur at excep- 
tional points of square root branching type. For certain parameter configurations two or 
more such exceptional points may happen to coalesce and to form a higher order branch 
point. We study the coalescence of two square root branch points semi-analytically for a 
VT— symmetric 4x4 matrix toy model and illustrate numerically its occurrence in the 
spectrum of the 2x2 operator matrix of the magneto- hydrodynamic a 2 — dynamo and of 
an extended version of the hydrodynamic Squire equation. 
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1 Krein space related physical setups and spectral phase transitions 

Some basic spectral properties of the Hamiltonians of VT— symmetric Quantum 
Mechanics (PTSQM) pQ can be easily explained from the fact [3] that these 
Hamiltonians are self-adjoint operators in Krein spaces |S] — Hilbert spaces 
with an indefinite metric structure. In contrast to the purely real spectra of self- 
adjoint operators in "usual" Hilbert spaces (with positive definite metric structure), 
the spectrum of self-adjoint operators in Krein spaces splits into real sectors and 
sectors with pair- wise complex conjugate eigenvalues. In physical terms these two 
types of sectors are equivalent to phases of exact VT— symmetry and spontaneously 
broken VT— symmetry pp. The reality of the spectrum of a PTSQM Hamiltonian, 
e.g., with complex potential ix 3 , means that this spectrum is located solely in a 
real sector and that the operator is quasi-Hermitian 1 ) in the sense of Ref. [7]- In 
general, the spectrum of a self-adjoint operator in a Krein space is spreading over 
both sectors. 

*) Presented at the 3rd International Workshop "Pseudo-Hermitian Hamiltonians in Quantum 
Physics" , Istanbul, Turkey, June 20-22, 2005. 

1 ) Because quasi-Hermitian operators form only a restricted subclass of self-adjoint operators 
in Krein spaces (pseudo-Hermitian operators in the sense of Ref. the question for the re- 

ality of the spectrum seems up to now only partially solved. Apart from the requirement for 
existing VT— symmetry of the differential expression of the operator, a subtle interplay between 
the operator domain and some, in general not yet sufficiently clearly identified, additional struc- 
tural aspects of its differential expression seems to be responsible for the quasi-Hermiticity of the 
operator (exact VT— symmetry of the corresponding PTSQM setup). 
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Considering the spectrum as a Riemann surface depending on model param- 
eters (not necessarily moduli) from a space M, the boundaries between real and 
complex spectral sectors can be described as algebraic variety ((multi-component) 
hypersurface) T in M. , T C M , where the Riemann surface has exceptional points 
of branching type [SJ E3 HH H2I 1131 HI) . The most generic varieties of this type 
are those of codimension one, codim(T) = dim(A4) — dim(T) = 1, what can be 
easily read off, e.g., from a pseudo-Hermitian 2 x 2— matrix model [151 1161 IT71 118) 

where T is simply a double cone y 2 = \w\ 2 C M 3 (y, 3tw, 3w). Such codimension- 
one varieties correspond to a pairwise real-to-complex transition of two eigenvalues 
what in Riemann surface terms is equivalent to a square root branching of two of its 
sheets (two spectral branches). At the same time the two (geometric) eigenvectors 
coalesce into a single (geometric) eigenvector, and an additional associated vector 
(algebraic eigenvector) appears. Instead of two eigenvalues of geometric and alge- 
braic multiplicity one (diagonal block) , a single eigenvalue of geometric multiplicity 
one, (m(A) = 1), and algebraic multiplicity two (n(A) = 2), (2 x 2— Jordan block) 
forms 18 . When the parameters on this codimension-one variety are further tuned 
to y = 8tw = 3ui = JS) one arrives at a diabolic point GO) where the single 
eigenvalue gets geometric and algebraic multiplicity two (m(X) = n(X) = 2, diag- 
onal 2 x 2— block). In a (more) general setting such configurations correspond to 
codimension-three varieties [181 l2T)j . 

Apart from these generic real-to-complex transitions on codimension-one vari- 
eties, there may occur higher order intersections of more than two Riemann sheets 
simultaneously — on varieties T of higher codimension, codim(T) > 2, and with 
larger Jordan blocks in the spectral decomposition [5J 1211 |2"2") . 

Below we present explicit examples for third order intersections in three Krein 
space related physical models. First, we sketch a few aspects of such intersec- 
tions semi-analytically (algebraically) for a maximally simplified VT~ symmetric 
4 x 4— matrix toy model. Afterwards, we show numerically that such intersections 
also occur in the operator spectra of the spherically symmetric a 2 — dynamo of mag- 
netohydrodynamics (MHD) [231 and of the recently analyzed VT— symmetric 
interpolation model of Ref. |2H) (which can be understood as a PT— symmetrically 
extended, rescaled and Wick-rotated version of the Squire equation of hydrody- 
namics) . 

2 Spectral triple points in a pseudo-Hermitian 4 x 4— matrix toy model 

Similar to the two different multiplicity contents (Jordan structures) of two- 
fold degenerate eigenvalues (see, e.g., also |H>J), one has to distinguish the following 
three types of Jordan structures/geometric multiplicities m(A) corresponding to a 
spectral triple point: 

— type I: m(A) = 1, n(A) = 3; 3 x 3— Jordan block; e.g., coalescence of two 
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square-root branch points connected by a purely real spectral segment (see 
Eq. (PJ and Figs. Q- El below), 

— type II: 771(A) = 2, 77(A) = 3; one 2 x 2— Jordan block + one simple eigen- 
value; e.g., a square-root branch point accidentally coincides with a simple 
eigenvalue, 

— type III: 777(A) = 3, 77(A) = 3; diagonal matrix A/3; e.g., three accidentally 
coinciding eigenvalues (generalized diabolic point). 

Subsequently, we concentrate on the physically most interesting case of type I triple 
points which correspond to pair-wise coalescing square-root branch points. 

As maximally simplified toy model we choose a 'PT--symmetric (pseudo-Her- 
mitian) 4 x 4— matrix setup in a diagonal representation of the parity (involution) 
operator V 

" = «^ = U-1 "«--)' P =(o -l)< (2 > 

where the 2 x 2— blocks satisfy the conditions H±± = H± ± , H^ = —H_, . In 

order to keep the demonstration as simple as possible, it suffices to consider a model 
based on two real (truncated) elementary pseudo-Hermitian 2 x 2— matrices 

tt _( 2=1,2 +2/1,2 ^1,2 \ _ „ ,„ c1B /on 

^1,2— ) £1,2,2/1,2, Wl o S K (6) 

as subsystems, which we embed VI — symmetrically into the also highly reduced 2 ) 
and real pseudo-Hermitian 4 x 4— matrix H with block structure © 



Hi , H-2 H — 



( xi+yi wi z \ 

x 2 + U2 w 2 

—wi x 1 — 7/1 

\ -z -w 2 x 2 - 7/2 / 



(4) 



The interaction between the subsystems Hi and H 2 is controlled by the coupling 
parameter z. For vanishing coupling, z = 0, the characteristic equation A(A) = 
det(if - A/ 4 ) = factors as A(A) = det(i?i - A/ 2 ) dct^ - A/ 2 ) = 0, and the four 
eigenvalues of the matrix H are defined by the eigenvalues of Hi and H 2 



\i,±=x 1 ±yy'(-wl, A 2:± = x 2 ± \Jyl - w\ . (5) 

The corresponding two square root branch points Ai = Xi, A2 = x 2 with A(A) = 
0, <9\A(A) = are located on the two (reduced) double cones 3 ) (crossed lines) 
defined by 

2 2 2 2 ia\ 

yj_=wi, y 2 = w 2- (6) 

2 ) Nine of the 16 effective free real parameters of the pseudo-Hermitian 4 X 4— matrix H are set 
to zero. 

s ) The two diabolic points 19 20 (with eigenvalues of geometric and algebraic multiplicity two) 
of the subsystems Hi , H2 are located at y\ = W\ = and j/2 = 7"2 = 0, respectively. 
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(7) 



For non- vanishing interaction z / 0, the eigenvalues of H are given as roots of the 
quartic equation 

A(A) = A 4 + a 3 A 3 + a 2 A 2 + ai A + a Q = 

with coefficients 4 ) 

a 3 = -2(xi+x 2 ) 

a 2 = -{y\ - w\) - (y\ - w\) + {x\ + x 2 ) 2 + 2x\x 2 + z 2 
oi = 2 [xi(y 2 - w\ - x\) + x 2 {y\ - w\ - x\)\ - z 2 (xi - y\ + x 2 + y 2 ) 



«o 



w 2 



x 2 ) + z 2 (xi - yi){x 2 + 3/2) 



(8) 



A smooth branch point coalescence can be easily arranged by VT— symmetrically 
blowing up a triple root 



Ai — A2 



(9) 



A3 — : /3(c) 7^ ^4( c ) 

("inflection point" configuration A(A) = <9aA(A) = <9 2 A(A) = 0) of the quartic 
equation J2J. For our illustration purpose, it suffices to explicitly derive a suit- 
able parametrization for an arbitrarily chosen triple root 5 ) satisfying the additional 
branch point conditions for the subsystems. Using these conditions together 
with ©, © in Newton's identities (Vieta's formulas) |2H1 

Ai + A 2 + A 3 + A 4 = — a 3 
AiA 2 + A1A3 + A1A4 + A2A3 + A2A4 + A3A4 = a 2 

A1A2A3 + A2A3A4 + A1A2A4 + A1A3A4 = —ai 

A1A2A3A4 = a (10) 



and setting for convenience /3( c ) = Z( c ) = 1, xi( c ) 
maining root A4( c ) and the three parameters X2(c)j 2/i,2(c) a s 



we easily obtain the re- 



A4(c)e 
^2(c)e 



3 + 2 3/2 e>0, e:=±l, 
^(A 4( c)c + 3) = 3 + 2 1 / 2 e, 

(3A 4 ( c ) e + 1 



<W9A 2 (c)£ 



2A 



4(c)e 



±1. 



A 



4(c)e 



f 



2V 9A 4(c) £ 



2A 



4(c)e 



1. 



(11) 



A computer algebraic test shows that the Jordan normal form of H for this triple 
root configuration reads 



H(c)j = 



I f f \ 

Oil 

1 

V A 4(c)e ) 



(12) 



4 ) The reason for restricting to the maximally simplified (truncated) 4 X 4— matrix setup TO 
was in keeping these coefficients aj., k = 0, . . . ,3 in a form sufficiently simple for quick inspection. 

5 ) A detailed discriminant-based 1271 case analysis and complete classification of possible inter- 
section scenarios for the four roots A fe , k = 1, , 4 will be presented elsewhere. 
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Fig. 1. The four solutions Ai,...,A4 of the quartic equation A(A) = (purely real 
branches highlighted fat) for different subsystem couplings z. For increasing z two of the 
exceptional points coalesce and disappear [at the inflection point of the inverted purely 
real branch t(5RAi): 9^ Al t(5RAi) = 0]. Simultaneously, the originally existing real spectral 
segment between the exceptional points disappears and the two complex valued sectors 

merge into a single one. 



so that the eigenvalue 0r c \ — 1 is indeed a type I triple point with geometric multi- 
plicity one and algebraic multiplicity three. For the 2D— illustration of the coalesc- 
ing branch points in Fig. ^ we have chosen a one-parameter matrix parametrization 
of the type H(t) = Hm + h(t), h(t — > 0) — > providing the blowing-up of the 
triple root for t ^ with 

xi = a:i( c ) , x 2 = x 2 { c ) + t, 

Vi = Vi(c)++ + 2i 2 ?/2 = 2/2(c)++ - t 3 , 

w\ = ?A( C )++ ~t w 2 = V2(c)++ + 3t 2 . (13) 

A characteristic feature of the present coalescing branch point setup consists in the 
merging of two originally separated complex spectral sectors (of broken VT— sym- 
metry) into a larger single sector and the simultaneous disappearance of the real 
segment between the two square-root branch points. This is in obvious contrast to 
coalescing branch point setups in simpler pseudo-Hermitian 2 x 2— matrix models, 
where the branch points are connected by two real or complex branches which 
disappear upon branch point coalescence. 



Czech. J. Phys. 55 (2005) 



5 



Uwe Giinther et al. 



3 Third order spectral transitions of the a 2 dynamo and the 

extended Squire setup 

Based on the information of the previous section, it is easy to identify con- 
figurations with coalescing branch points in the spectra of the MHD a 2 — dynamo 
operator and of the operator of the extended Squire setup. 

In case of the operator [181 1241 l2"B*] of the spherically symmetric a 2 — dynamo 

m 

HM = 



-Q[i] 



Q[c] := -(d r + l/r)a(r)(d r + l/r) + a (r)^4^ 



Q[a] -Q[l] 

with a— profile 6 ) 

a(r) = C [-(21.465 + 2.4670 + (426.412 + 167.928C)r 2 
-(806.729 + 436.289C)r 3 + (392.276 + 272.991C)r 



(14) 



(15) 



and physical (realistic) boundary conditions [181 1231 |2~H 125) a triple point transition 
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Fig. 2. Coalescing branch points in the spectrum of an a 2 — dynamo with a— profile 1 1 U 
depending on a warp parameter £. 

in the decay-mode sector (5RA < 0) is depicted in Fig. [3 Up to now it is still an 
open question whether such a transition in the sector of growing modes (5RA > 0) 
will have any physical significance, e.g., in dynamo experiments |3(JI 131) . 

A different situation occurs in the case of the extended Squire setup of Ref. 

[-d 2 y +gy 2 {ivY} ^(y) = E ^(y), ^(y = ±b) = (16) 

6 ) The rather special numerical coefficients in the a— profile are adopted from the recently 
studied field-reversal scenario for a 2 — dynamos 1291 . 
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which corresponds to the Bender-Boettcher problem pQ over a finite interval [—6, 6]. 




Fig. 3. Coalescing branch points in the spectrum of the extended Squire setup 1161 with 
unit coupling, g = 1, in dependence of the cut-off /box length b. 

There a coalescence of two square-root branch points occurs close to the real-to- 
complex transition point in the associated Herbst-box model [23 (Herbst model 
over a finite interval) located at the low-energy end of a web-like branch structure. 
An example for the corresponding triple point transition (which supplements the 
qualitative and graphical analyzes of Ref. [35]) is depicted in Fig. 01 

This work has been supported by DFG grant GE 682/12-2. 
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